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ABSTRACT 

We propose a model of magnetic connection (MC) of a black hole with its sm-rounding 
accretion disc based on large-scale magnetic field. The MC gives rise to transport of 
energy and angular momentum between the black hole and the disc, and the closed field 
lines pipe the hot matter evaporated from the disc, and shape it in the corona above the 
disc to form a magnetically induced disc-corona system, in which the corona has the 
same configuration as the large-scale magnetic field. We numerically solve the dynamic 
equations in the context of the Kerr metric, in which the large-scale magnetic field is 
determined by dynamo process and equipartition between magnetic pressure and gas 
pressure. Thus we can obtain a global solution rather than assuming the distribution 
of large-scale magnetic field beforehand. The main MC effects lie in three aspects. (1) 
The rotational energy of a fast-spinning black hole can be extracted, enhancing the 
dissipation in the accretion disc, (2) the closed field lines provide a natural channel 
for corona matter escaping from disc and finally falling into black hole, and (3) the 
scope of the corona can be bounded by the conservation of magnetic flux. We simulate 
the high-energy spectra of this system by using Monte-Carlo method, and find that 
the relative hardness of the spectra decreases as accretion rate or black hole spin 
increases. We fit the typical X-ray spectra of three black-hole binaries (GRO J1655-40, 
XTE 1118+480 and GX 339-4) in the low/hard or very high state. 

Key words: accretion, accretion discs — black hole physics — X-rays: binaries — 
magnetic fields 



1 INTRODUCTION 

It is well known that black hole binary transients present 
different spectra in different stages during their outbursts, 
which are usually referred to as low/hard state, high/soft 
state and intermediate (or very high) state (McClintock & 
Remillard 2006, hereafter MR06). It is commonly realized 
that how to interpret the different spectral profiles and the 
inducement of state transitions may be the key to look into 
the real physical processes during the outbursts. 

Accretion disc is generally considered as the main en- 
ergy source for compact objects. Historically, people pro- 
posed various accretion models to interpret different ob- 
served spectra and light curves, such as standard thin disc 
(Shakura & Sunyaev 1973, hereafter SS73; Pringle 1981) and 
advection dominated accretion flow (ADAF, Narayan & Yi 
1994). The main criterions to define different spectral states 
are the relative weight of power-law component to ther- 
mal component in X-ray spectra, and the photon index of 
power- law component (MR06). However, we cannot expect 
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to interpret the observed various spectra based on a single 
accretion mode. For example, standard thin disc is largely 
good for soft spectra of high/soft state, but it cannot in- 
terpret the hackneyed power-law component in high energy 
spectra. ADAF is congenitally propitious to hard spectra of 
low/hard state, but it is invalid in fitting the observed soft 
excess and also inverse radio spectra. Recently, it is generally 
considered that different spectral components have different 
origins, i.e., the thermal component of X-ray spectra comes 
from the black-body radiation of a standard thin disc, while 
the power-law component arises from corona or jet. And ra- 
dio emission is usually considered to be dominated by the 
synchronization radiation from jet. 

Angular momentum transport is the essence of accre- 
tion theory. Shakura & Sunyaev (1973) proposed the semi- 
nal alpha viscosity law to discuss complicated viscous pro- 
cess. Balbus & Hawley (1991) argued that magnetic rotation 
instability of tangled magnetic field within accretion disc 
would be the physical origin of viscosity, and it is generally 
regarded as another milestone of accretion theory. On the 
other hand, the heating mechanism for corona is another 
open question, in which magnetic field plays an essential 
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role. The most popular sketch is that the tangled small-scale 
magnetic field amplified by dynamo process buoys up into 
corona and then reconnects, heating corona by magnetic en- 
ergy dissipation. Meanwhile, the collimation mechanism of 
observed rclativistic jots is also under controversy. Recently, 
Miller et al. (2008) carefully investigated the origin of disc 
wind of GRO J1655-40 during its outburst in 2005. Inter- 
estingly, he ruled out several popular sketches for launching 
disc wind and concluded that the disc wind of GRO J1655-40 
must bo driven by magnetic field. More and more attention 
has been paid to the important role of magnetic fields in 
spectral states of black hole binaries. 

As is well known, the Blandford-Znajeck (BZ) process is 
an important mechanism for jet production, in which energy 
is extracted efficiently from a spinning black hole to power 
remote astrophysical loads by invoking "open" large-scale 
magnetic field and frame dragging effect. The BZ process can 
be regarded as the magnetic Penrose process (Penrose 1969; 
Blandford & Znajeck 1977; Livio, Ogilvie & Pringle 1999). 
As a variation of the BZ process, closed field lines connecting 
a black hole with its surrounding accretion disc lead to the 
transfer of energy and angular momentum between the hole 
and the disc, and it is referred to as the Magnetic Connection 
(MC) process (Li 2002; Wang et al. 2002). The direction of 
such transport is determined by the angular velocity of the 
hole relative to that of the disc. Energy can be extracted 
efficiently from a spinning black hole with a* > 0.3594 to 
the inner accretion disc in the MC process, enhancing disc 
radiation and a very steep emissivity index (Wilms et al. 
2001; Wang et al. 2003). 

However, the importance of the MC remains contro- 
versial. Tomimatsu & Takahashi (2001) obtained a global 
magnetospheric structure with a closed-loop and open field 
lines threading the inner and outer parts of the disc, and 
the poloidal magnetic field is generated by a toroidal elec- 
tric current in a thin disc with an inner edge based on 
the vacuum Maxwell equations in the Schwarzscliild back- 
ground. They found that the MC between the black hole 
and the disc breaks down if a uniform external field is strong 
enough. Hirose et al. (2004, hereafter H04) presented a de- 
tailed analysis of the magnetic field structure based on a set 
of three-dimensional general rclativistic magnetohydrody- 
namic (GRMHD) numerical simulations of accreting tori in 
the Kerr metric with different black hole spins. Very recently. 
Fragile & Meier (2008, hereafter FM08) presented one of 
the first physically-motivated two-dimensional GRMHD nu- 
merical simulations of a radiatively-cooled black hole accre- 
tion disc, and they obtained a magnetically-dominated ac- 
cretion flow (MDAF) in the interior with an outer advection- 
dominated accretion flow (ADAF). It turns out that the im- 
portance of the MC depends on the radiative cooling, which 
is not taken into account in the non-radiative simulations 
of H04. As argued in FM08, a more ordered magnetic field 
could be formed due to the radiative cooling, by which a 
dramatic increase in the dominance of magnetic stresses is 
created. Considering the fact that a standard thin disc with 
a very effective radiative cooling is adopted in our model, 
we expect that the MC configuration of the closed magnetic 
field lines can be formed. 

Another controversial issue involves the stresses on the 
inner edge of an accretion disc. Some authors (Krolik 1999; 
Gammie 1999; Agol & Krolik, 2000) argued that the torques 



exerted on the inner edge of a disc might arise from the MC 
to the plunge region. Hawley and Krolik (2002) showed that 
magnetic stresses between the plunge region and the inner 
disk grow 10 times the value of stresses in the disk body. 
Garofalo & Reynolds (2005) showed the importance of such 
inner disk torques and how they modify the dissipation pro- 
file of a standard accretion disk even for Schwarzschild black 
holes. On the other hand, Paczynski (2000) and Afshordi & 
Paczynski (2003) suggested that the zero-torque condition 
is likely to be a good approximation for geometrically thin 
discs. This result was confirmed by Shafee et al. (2008) who, 
using a global height-integrated model, showed that modi- 
fications to the stress profile are negligibly small for disk 
thicknesses h less than about a tenth of the local radius r. 
Since only a thin disc with corona is involved in our model, 
we neglect the stresses on the inner edge of the accretion 
disc for simplicity. 

Corona is a layer of hot tenuous plasma with temper- 
ature Te >~ lif K above accretion disc. The soft photons 
emitted from disc surface are scattered up by relativistic 
electrons in corona, resulting in a much harder spectrum. 
Corona is usually invoked to interpret the power-law 
spectral component from tens of KeV to hundreds of 
KeV observed in black hole binaries. The main differences 
between various corona models lie in the heating mechanism 
and the configuration of corona (i.e., temperature, density 
and geometry). Stella & Nosner (1984) argued that tangled 
magnetic field amplified by dynamo process could float into 
corona and reconnect there, thus corona is heated by the 
dissipated magnetic energy. With the assumed corona ge- 
ometry, Liu, Mineshige & Shibata (2002, hereafter LMS02) 
proposed a disc-corona system by resolving equation of 
energy with the balance between thermal conduction and 
corona evaporation. It is noted that coronae make spectra 
harder, but cannot augment total luminosities. In addition, 
Merloni & Fabian (2002, hereafter MF02) pointed out that 
corona is also a perfect launching site for outflow from 
accretion disc. 

In this paper, we propose a magnetic disc-corona model, 
in which the magnetic field configuration is determined self- 
consistently, and we discuss the effects of magnetically in- 
duced corona on spectra of a relativistic thin Keplerian disc. 
Based on a few assumptions, we obtain a global solution of 
dynamic equations of accretion disc, by which the interac- 
tion of corona, the origin of large-scale magnetic fleld and 
the MC process are clarified. In this model, the MC pro- 
cess not only transfers rotational energy from black hole to 
the inner disc, but also provides a natural tube for evapo- 
rated hot plasma to fall into the hole as argued in FM08. 
The tenuous hot plasma is constrained and forced to move 
along the magnetic field lines, giving rise to corona natu- 
rally according to the configuration of the closed field lines. 
We assume that corona is of the same configuration as the 
closed magnetic field for the MC process, and calculate the 
boundary radius of corona based on evolving the magnetic 
field configuration. It turns out that the disc-corona sys- 
tem overwhelmingly dominates the high-energy radiation of 
whole system. 

Both the tangled small-scale magnetic field within ac- 
cretion disc and the ordered large-scale magnetic field re- 
lated to the MC process are involved in this model. It is as- 
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sumed that the energy density of small-scale magnetic field 
and the gas pressure within the disc are of the same order 
based on energy equipartition. Magnetic rotational insta- 
bility (MRI) or turbulence of tangled small-scale magnetic 
field is expected to work for viscous dissipation, while ex- 
traction and transport of energy and angular momentum 
arise mainly from the large-scale magnetic field. The small- 
scale magnetic field can be amplified by dynamo process, 
and corona is heated by magnetic buoyancy and reconnec- 
tion. Large-scale magnetic field is generated from the small- 
scale magnetic field, and the former is scaled by the disc 
radius and its intensity is usually 1 ~ 3 order weaker than 
the latter. We present a numeric code to solve the coupled 
dynamic equations of this magnetic disc-corona system self- 
consistently and simulate its spectra by using Monte-Carlo 
method. 

This paper is organized as follows. We outline the model 
in §2, where the coupled disc-corona system and the MC 
process are included. In §3, we numerically resolve the mag- 
netic disc-corona system and simulate the high-energy spec- 
tra from the inner disc by using Monte-Carlo method. And 
four sampling results are presented with detailed analysis. 
In §4, we fit the typical X-ray spectra of three black-hole 
binaries during their outbursts. §5 is the discussion with a 
brief conclusion. 

The Boyer-Lindquist coordinates and geometric units 
G = c = 1 are used throughout this paper. 



where Bd, Va and p are the intensity of interior tangled 
small-scale magnetic field, local Alfven speed and mass den- 
sity, respectively. Te// is the effective temperature of accre- 
tion disc and a is the Stefan-Boltzman constant. 

In the context of relativistic hot corona, inverse Comp- 
ton scattering is very effective, by which corona cooling is 
dominated. According to energy balance given by LMS02, 
Qtor = Qcomp, we have 



Bi 



(4) 



where Urad ~ '^T^ff ~ ^Qd radiant energy density at 

vicinity of disc surface. Tcor and Tcor are the temperature 
and the optical depth of corona. The quantities k, irie and a 
are the Boltzman constant, electron mass and radiation con- 
stant, respectively. Observations reveal that optical depths 
of corona for AGNs and for black hole binaries in low/hard 
state (LHS) lie in a very narrow range around Tcor ~ 1 (Gier- 
linski et al. 1997; Zdziarski 1999), and we adopt Tcor ~ 1 for 
facility, i.e., 



Tcor = UcorTTl ~ 1, 



(5) 



where n^or and / are the number density of electron and the 
height of corona, respectively. 

2.2 The MC Process 



2 MODEL DESCRIPTION 
2.1 Disc and Corona 

In standard thin disc model, the interior viscous stress tr^, is 
usually assumed to be proportional to total pressure Ptot (in- 
cluding gas pressure, radiation pressure and magnetic pres- 
sure), namely the famous alpha viscosity law (SS73). Vis- 
cous stress is alternatively assumed to be proportional to gas 
pressure Pgas, or radiation pressure Prad or the geometrical 
mean of gas pressure and total pressure ^ Pgas ■ Ptot and so 
on (Stella & Nosner 1984; Wandel & Liang 1991; Taam & 
Lin 1984; Wang et al. 2004, table 2 and references therein). 
Wang et al. (2004) fitted the X-ray luminosities of 56 radio- 
quiet AGNs and concluded that the observations prefer to 
the interior viscous stress proportional to gas pressure. If 
interior viscous process is dominated by tangled small-scale 
magnetic field, the viscous pressure should be comparable 
to magnetic pressure (Balbus & Hawley 1991; MF02), and 
we assume 



OiPgas 
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(1) 



According to typical disc-corona scenario, part of the 
viscously dissipated energy Q is released as in the disc, 
emitting eventually as black-body radiation and supplying 
seed photons for Comptonization of corona. The rest dis- 
sipated energy, Qtor, heats corona and maintains its rela- 
tivistic temperature via magnetic reconnection. The quan- 
tity Qtor is proportional to magnetic energy density and 
local Aflven speed (MF02; LMS02), and we have 



Q — Qd ~i~ Qcor — Q cor I f cor 



) + 

C cor 



Bf 



(2) 
(3) 



Being different from standard thin disc, energy dissipation 
Q in this model arises from two sources: one is the gravita- 
tionally bound energy of accreting matter, and the other is 
the rotational energy extracted from the black hole via the 
MC process (Wang et al. 2002). 

A fast-spinning black hole exerts a significant torque on 
its surrounding disc, if the two are coupled by the closed 
magnetic field. Energy and angular momentum are trans- 
ferred between the two simultaneously. The torque exerted 
by the black hole onto an infinitesimal annular region of the 
disc is given by 



dTMC = 2 



(f) 



dZH 



= Airrdr ■ H 



MC, 



(6) 



where Hmc = (l/47rr-) ■ dTMc/dr is the fiux of angular 
momentum transferred from the black hole to the disc. The 
quantities Qh and Qd are the angular velocity of the black 
hole and that of the disc, respectively. The quantity tp = 
tp{r, 6) is the magnetic fiux through a surface bounded by a 
curve with r = constant and 6 = constant, and we have 



dtp = Bp-2-Kr 




=^/2= -Bh-2'KTh-2M sinOde |r=r„, (7) 



th is the radius of black-hole horizon and the resistance of 
the horizon corresponding to the magnetic flux is written as 
(Wang et al. 2002) 



dZ, 



(Rh\ th + M'^al cos^ 



V 27r y {r]j + NPal)8ine 



(S 



• dr 



(8) 



Rh is the surface resistivity of black-hole horizon. M and 
a, are the mass and the dimensionless spin of black hole, 
respectively. 

In equation ([7]), Bp is the poloidal component of the 
large-scale magnetic field anchored on the disc. Equation ((Tj 
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is derived from conservation of magnetic flux, by which the 
mapping relation r{9) between black-hole horizon and the 
inner region of the disc can use determined. Assuming that 
the inner edge of the disc connects the horizon a,t = 7r/2, 
we obtain the outer boundary of the MC region, r{0 = 0), 
by integrating equation ((7} from 9 = n /2 to 0. 

Still, we have to estimate magnetic field strength on 
the horizon required by equations (|6} and ([7]). Assuming 
a uniform magnetic field on the horizon and adopting the 
equilibrium between magnetic pressure and the ram pressure 
given by Moderski, Sikora & Lasota (1997), we have 



Bh = j2Mn/rl, 



(9) 



where Mo is accretion rate. 

The relativistic conservation equations of energy and 
angular momentum for a disc with magnetic connection can 
be written as (Li 2002) 

d_ 

dr ^ 
d_ 
dr ^ 

The quantities Q and g are the dissipated energy per 
unit disc surface and interior viscous torque of the disc, re- 
spectively. And E\ and Qd are respectively the specific 
energy, specific angular momentum and angular velocity of 
the test particles in the equatorial plane of a Kerr black hole 
(Bardeen, Press & Teukolsky 1972). 



-{MdL'' - <?) = 47rr(QLt - Hmc), 
-{MdE'' - <? ■ nc) = 47rr(OSt - Hmc^d) 



(10) 
(11) 



no = ^M/r^ /(l + a,/f2), 



(l-2/f + a./f 2 ) 



(l-2al/f^+al/f^) (12) 



•\/ \-3/f + 2a,/f 2 

where f = r/M. The quantities A and "D in equation ((Tj are 
relativistic correction factors in Kerr metric, and they read 



A = l + ajr + 2a: /r 



V 



2/f + al/r 



(13) 



We can derive the energy dissipation Q and interior viscous 
torque g by resolving equations pO|l and (|lip . 



Q{r) 



{da + Qmc = Qda+ 



7 (^"^[)2 ■ IL - i^nL^)HMcrdr 



(14) 



Qda ~ 



1 



—dQo/dr 



4nr (El - floL^y 



f {e'' - !^DLt)Afo-dLt(i5) 



-4:7vr ■ Q{r) 



(16) 



-dQo/dr 

rms is the radius of the innermost stable circular orbit, which 
is assumed to be the inner edge of accretion disc. 

In this model we intend to resolve the dynamic equa- 
tions combining the MC effects, which might change the 
dynamic property of accretion disc significantly. To resolve 
equations ([T4|) - (fT6|) . we adopt the relation between the large- 
scale magnetic field above the disc and tangled small-scale 
magnetic field in the disc as follows (Livio et al. 1999), 



(17) 



where h is the half height of the disc. The intensity and 
distribution oi Bd and Bp can be derived by resolving the 
dynamic equations. 



3 NUMERICAL RESULTS 

For facility, we consider a slab corona with I ~ lOrma. 
Inspecting equations (|5))- p6p we find that magnetic con- 
nection at radius r could affect the disc region beyond r. 
We design an iterative algorithm and integrate equations 
(|10l) - (|16| ) from the innermost stable circular orbit (ISCO) 
to the outer boundary of the MC region, and obtain self- 
consistently global solutions of the disc-corona system, in 
which the full solutions of corona and the MC effects are 
included. As argued above, energy can be extracted from a 
fast-spinning black hole via the MC process, resulting in the 
enhancement of dissipation in the disc. Subsequently, the 
disc height, pressure and magnetic field increase. It is noted 
that the MC effects scaled by the square of magnetic field 
become very strong for the amplified magnetic field. Since 
corona heating via magnetic reconnection is scaled by the 
cube of magnetic field [cf. equation([3ll], we expect that a 
stable saturated magnetic configuration can be maintained 
in our model. 

The iterative algorithm consists of the following steps: 
(1) resolving numerically the disc-corona system without the 
MC process; (2) deriving the large-scale magnetic field and 
calculating the energy extraction by the MC process; (3) 
combining the energy dissipation of the disc with the MC 
process, and resolving the disc-corona system; (4) repeating 
steps (2) and (3) until the magnetic field is saturated; (5) 
integrating outward to a larger radius (r ~ r -f dr); (6) 
repeating (l)-(5) until the outer boundary of MC process ( 
6{r) = ) is reached. 

Here we present four typical solutions as shown in 
Figure 1-2, 3-4, 5-6 and 7-8, where Td is the temperature 
of accretion disc on the equatorial plane. The concerned 
parameters are listed in Table 1, in which the black-hole 
mass mth (= M/Mq), distance D, inclination angle i, 
corona height I, viscous coefficient a, black-hole spin a* 
and accretion rate m (in unit of the Eddington rate) are 
included. The first three parameters are constrained by 
observations. The viscous coefficient is taken as a typical 
value, a = 0.3. In fact the corona height might vary with 
the disc radius, and its value could influence the density 
of corona. Fortunately, emerged spectrum from corona is 
mainly determined by the temperature and optical depth. 
Since we adopt a constant optical depth for corona, our 
results are not sensitive to and we fix Z = 10r„is in 
calculations. So there are actually only two free parameters 
in our model, accretion rate rh and black-hole spin a*. 
The former represents the "fuel supply", and the latter 
determines the radiant efficiency (luminosity per unit 
mass). Both parameters are very important for the emerged 
spectrum. Actually, we can demonstrate that they are the 
key elements of the spectral profile. 

The MC effects on accretion disc. The MC process 
utilizes the difi^erential rotation between a black hole and its 
accretion disc to transfer energy and angular momentum. As 
argued by Wang et al. (2002) , the MC effects become much 
stronger for greater black hole spin. In this paper, we find 
that the MC effects on accretion disc also depend on accre- 
tion rate. The distribution of closed field lines and corona is 
very concentrated to the innermost region of accretion disc 
(<~ 20M) for great accretion rate (>~ 0.2MEdd)- Thus the 
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Table 1. Concerned Paraniotors in the samyjlin,t^ results and simulations 



Source 




Observation 






Model 




Results 






D /kpc i 


1 degree 


^ / "^ras 


OL 


a* 


m 


(Figure) 


Sample 1 


10.0 


3.0 


0.0 


10.0 


0.3 


0.900 


0.300 


1, 2 


Sample 2 


10.0 


3.0 


0.0 


10.0 


0.3 


0.200 


0.030 


3, 4 


Sample 3 


10.0 


3.0 


0.0 


10.0 


0.3 


0.200 


0.300 


5, 6 


Sample 4 


10.0 


3.0 


0.0 


10.0 


0.3 


0.900 


0.030 


7, 8 


GRO J 1655-40 (LHS) 


7.0" 


3.2'' 


70.0= 


10.0 


0.3 


0.200 


0.030 


9 


GX 339-4 (VHS) 


5.8'' 


7.0" 


12.0-'' 


10.0 


0.3 


0.600 


0.250 


10 


XTE J1118-I-480 (LHS) 


8.09 


1.8'' 


70.0' 


10.0 


0.3 


0.400 


0.005 


11 
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et al., 1998, MNRAS, 301, 435; (f) Miller J. M., Fabian A. C., et al., 2004, ApJ, 606, L131; (g, h) McClintock J. E., Garcia M. R., et 
al., 2001, ApJ, 551, 147; (i) Wagner R. Mark, Foltz C. B., et al., 2001, ApJ, 556, 42. 



energy pumped into the disc via the MC process is compa- 
rable to the released gravitational energy and the radiation 
efficiency of the system is enhanced significantly. As shown 
in Figure la and lo, both temperature and dissipation rate 
of the disc are augmented obviously due to the MC effects. 
It is interesting to note that the MC process has little influ- 
ence on gas pressure, but it strengthens radiation pressure 
significantly. As a result, the inner region of disc becomes 
thick enough, being apt to form large-scale magnetic field. 
On the other hand, the stronger large-scale magnetic field is, 
the more efficient is MC process. So the energy dissipation 
and the MC process promote each other in our model. When 
accretion rate is small (<~ 0.05MBdd), the distribution of 
closed field lines and corona can extend to >~ lOOM. In 
this case, the magnetic field is too weaJj to give rise to ob- 
vious efi^ects on the disc dynamics, and the closed field lines 
mainly provide a channel for corona matter escaping from 
the disc and finally falling into the hole. 

Recalling our previous researches, we conclude that 
transport of energy and angular momentum via the MC 
process is generally concentrated in a very narrow region 
(~ lOM) near the inner disc boundary. However, it can 
change the total luminosity of the whole system significantly. 

Emerged spectra of the disc-corona system. Wc 

find that spectral profiles change obviously with both ac- 
cretion rate and black-hole spin. Generally, thermal radiar 
tion mainly comes from the disc, and power-law component 
comes from corona, implying that the higher corona tem- 
perature has the harder spectrum. 

(1) The fraction of thermal component increases as ac- 
cretion rate increases. When accretion rate increases, the 
magnetic field and thermal production are enhanced, result- 
ing in much more soft photons escaping from the surface of 
disc. Thus the corona temperature is suppressed effectively 
by the cooling of Comptonization. 

(2) The fraction of power-law component decreases as 
black hole spin increases. The gravitational potential well 
becomes very deep as the black hole spins fast. On the other 
hand, the stronger frame dragging effects give rise to more 
energy transferred from a spinning hole to the disc in the 
MC process. Both effects can strengthen the radiation from 
the inner disc as well as suppress the corona temperature. 

Generally speaking, the total luminosity of the disc- 



corona system is positively correlated to accretion rate and 
black-hole spin. We find that the spectra are always domi- 
nated by thermal component when accretion rate or black- 
hole spin is large (e.g., Md > O.SMedd, i* > 0.9), and these 
results are in agreement with those obtained by Esin, Mc- 
CUntock & Narayan (1997). 



4 COMPARISON TO OBSERVATIONS 

We try to fit the typical X-ray spectra of 3 black-hole bina- 
ries in the low/hard or very high state, i.e., GRO J1655-40, 
XTE 1118-1-480 and GX 339-4. The data are taken from lit- 
eratures (without error bars): GRO J1655-40 (Brocksopp et 
al. 2006), XTE 1118-^480 (Yuan et al. 2005 and references 
therein), GX 339-4 (Miller et al. 2004). 

As shown in Figures 9 and 11, the hard X-ray spectra 
with low luminosity (GRO J1655-40 and XTE 1118-1-480 in 
low/hard state) arc well fitted. In these simulations, the ac- 
cretion rates are small, gravitational binding energy is mildly 
released along the inner region of accretion disc and most 
of the energy (larger than 80%) is dissipated in corona via 
magnetic reconnection (cf. Figure 3f). And it turns out that 
the spectra are hard but dim. Inspecting Figure 10, we find 
that the soft X-ray spectrum with very high luminosity of 
GX 339-4 in very high state can also be well fitted. In this 
simulation, the accretion rate and black-hole spin arc large. 
Gravitational binding energy is released efficiently and effi- 
cacy of MC process also becomes significant. Energy dissi- 
pation is concentrated at the inner most region of the disc. 
Meanwhile, there is only a small fraction of dissipated en- 
ergy distributed to corona (cf. Figure If). And it turns out 
that the spectrum is soft but bright. 

However, there are excesses in soft X-ray band in these 
simulations, when we try to fit the total X-ray spectra. Since 
it is an embryo model, we cannot expect our model to fit the 
full features of the spectra. 



5 DISCUSSION 

In this paper, wo propose a magnetic disc-corona model, 
in which the configuration of large-scale magnetic field and 
corona can be self-consistently determined with several min- 
imal assumptions. We combine accretion disc, corona, origin 
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of large-scale magnetic field and the MC process into this 
model to investigate the interaction among them. The global 
solutions are obtained numerically. The size of corona can 
be self-consistently determined by resolving the magnetic 
configuration of the MC process. It is found that the MC 
process could change the properties of accretion disc signif- 
icantly for a large accretion rate. In addition, we simulate 
the emerged spectra from the inner region of accretion disc 
by using Monte-Carlo method. 

Observations reveal that large-scale magnetic fields ex- 
ist in compact objects. However, the origin of large-scale 
magnetic field is still under controversy. The popular treat- 
ment of large-scale magnetic field is invoked magnetohy- 
drodynamical (MHD) simulations. It is particularly noted 
that Uzdensky (2004, 2005) investigated the MC between a 
black hole and its surrounding disc, and he solved the Grad- 
Shafranov equations to determine the configuration of large- 
scale magnetic field in a force-free magnetosphere. Instead, 
we use the phenomenal description of dynamo mechanism 
to avoid complex calculations, and treat the origin of large- 
scale magnetic field directly by invoking the dynamics of 
accretion disc. In our model, the large-scale magnetic field, 
corona and accretion disc are connected very closely, which 
cannot be separated from each other. 

It turns out that in this model the emerged spectra 
change significantly with accretion rate, naturally giving 
rise to hard spectra with low luminosity and soft spectra 
with high luminosity. This feature might be related to the 
state transitions of black-hole binaries. As argued by Esin 
et al.(1997), the state transition could be induced by the 
variation of accretion rate. As a matter of fact, the key for 
state transitions is the accretion rate of the inner disc region 
rather than the rate of mass transferred from a donator star. 
The variation of the accretion rate of the inner disc region 
mainly has two possibilities: (1) mass rate from donator star 
changes; (2) a large fraction of accreting matter is thrown 
off the disc, e.g. by the Blandford-Payne (BP) process, in 
the form of jet or disc wind (Blandford & Payne 1982). It 
seems reasonable that most of accreting matter escapes from 
accretion disc before it drops into black hole in the low/hard 
state of black-hole binaries. 

The BP process is a very important mechanism for 
disc wind and jet production, which is also related directly 
to large-scale magnetic field. In the BP process, accreting 
matter could be thrown centrifugally off accretion disc 
and accelerated along the field lines. The outflow can be 
self-coUimated to jet due to magnetic pinch effect. However, 
it would become disc wind for weak magnetic field. The 
BP process can change the accretion rate of the inner 
disc region likely and thus affect the high energy emission 
of the system. It might play a very important role in 
state transitions of black-hole binaries. The outflow is also 
expected to produce some radiation and absorptive effects. 
Although there is little evidence of BP winds in GRMHD 
(see e.g. H04), the observed rich absorptive lines of GRO 
J1655-40 are likely related to winds driven by the BP 
process (Miller et al. 2008). The presence of BP process in 
our model would help to fit the observations and we leave 
it to the future work. 
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Figure 1. A global solution of the magnetic disc-corona system in the cases of large accretion rate and large black-hole spin (Sample 
1). Quantities are plotted in dashed and solid lines for disc -|- corona and disc + corona -|- MC, respectively. 
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Figure 2. Emerged spectrum from the disc-corona system in the case of large accretion rate and large black-hole spin (Sample 1). The 
total emissive spectrum and its thermal, comptonized and reflective components are plotted in thick-zigzag, solid, dashed and thin-zigzag 
lines, respectively. 
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Figure 3. A global solution of the magnetic disc-corona system in the cases of small accretion rate and small black-hole spin (Sample 
2). Quantities are plotted in dashed and solid lines for disc -|- corona and disc -|- corona -|- MC, respectively. (The dashed and solid lines 
almost coincide completely.) 
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Figure 4. Emerged spectrum from the disc-corona system in the cases of small accretion rate and small black-hole spin (Sample 2). The 
total emissive spectrum and its thermal, comptonized and reflective components are plotted in thick-zigzag, solid, dashed and thin-zigzag 
lines, respectively. 
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Figure 5. A global solution of the magnetic disc-corona system in the cases of large accretion rate and small black-hole spin (Sample 
3). Quantities are plotted in dashed and solid lines for disc -|- corona and disc -|- corona -|- MC, respectively. (The dashed and solid lines 
almost coincide completely.) 
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Figure 6. Emerged spectrum from the disc-corona system in the cases of large accretion rate and small black-hole spin (Sample 3). The 
total emissive spectrum and its thermal, comptonized and reflective components are plotted in thick-zigzag, solid, dashed and thin-zigzag 
lines, respectively. 
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Figure 7. A global solution of the magnetic disc-corona system in the cases of small accretion rate and large black-hole spin (Sample 
4). Quantities are plotted in dashed and solid lines for disc -t- corona and disc + corona + MC, respectively. 
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Figure 8. Emerged spectrum from the dise-corona system in the cases of small accretion rate and large black-hole spin (Sample 4). The 
total emissive spectrum and its thermal, comptonized and reflective components are plotted in thick-zigzag, solid, dashed and thin-zigzag 
lines, respectively. (Note: The comptonized and reflective components are too weak and absent from the flgure) 
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Figure 9. GRO J1655-40 in low/hard state (2005.03.06). The Concerned parameters are listed in Table 1. The embedded small figure 
is the emerged spectrum of our model. It is shown that the hard X-ray spectrum ( >~ 2 KeV) can be well fitted. It is expected that the 
excess in soft X-ray band could be absorbed by the surrounding medium (cf. Brocksopp et al. 2006). 
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Figure 10. GX 339-4 in very high state (2002.09.29). The Goncciiicd parameters are listed in Table 1. The embedded small figure is 
the emerged spectrum of our model. It is shown that the hard X-ray spectrum ( 2 KeV) can be well fitted. It is expected that the 
excess in soft X-ray band could be absorbed by the surrounding medium (cf. Miller et al. 2004). 
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Figure 11. XTE J1118+480 in low/haxd state. The Concerned parameters are listed in Table 1. The embedded small figure is the 
emerged spectrum of our model. It is shown that the X-ray spectrum ( >~ 0.1 KeV) can be fitted in the rough. The excess in lower 
energy band below 0.1 KeV might arise from bolometric radiation from the outer disc region and synchrotron radiation of jet, which is 
beyond the scope of this paper, (cf. Yuan et al. 2005) 



